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The synthesis and crystallographic and physical properties of poly crystalline EuNiGe3 are re- 
ported. EuNiGe3 crystallizes in the noncentrosymmetric body-centered tetragonal BaNiSn3-type 
structure (space group 1 4mm), in agreement with previous reports, with the Eu atoms at the cor- 
ners and body center of the unit cell. The physical property data consistently demonstrate that 
this is a metallic system in which Eu spins S = 7/2 order antiferromagnetically at a temperature 
Tn = 13.6 K. Magnetic susceptibility \ data for T > Tn indicate that the Eu atoms have spin 7/2 
with g = 2, that the Ni atoms are nonmagnetic, and that the dominant interactions between the Eu 
spins are ferromagnetic. Thus we propose that EuNiGe3 has a collinear A-type antiferromagnetic 
structure, with the Eu ordered moments in the ab-plane aligned ferromagnetically and with the mo- 
ments in adjacent planes along the c-axis aligned antiferromagnetically. A fit of x(T < Tn) by our 
molecular field theory is consistent with a collinear magnetic structure. Electrical resistivity p data 
from Tn to 350 K are fitted by the Bloch-Griineisen model for electron-phonon scattering, yielding a 
Debye temperature of 265(2) K. A strong decrease in p occurs below Tn due to loss of spin-disorder 
scattering. Heat capacity data at 25 K < T < 300 K are fitted by the Debye model, yielding the 
same Debye temperature 268(2) K as found from p(T). The extracted magnetic heat capacity is 
consistent with S = 7/2 and shows that significant short-range dynamical spin correlations occur 
above Tn. The magnetic entropy at Tn = 13.6 K is 83% of the expected asymptotic high-T value, 
with the remainder recovered by 30 K. 

PACS numbers: 75.50. Ee, 75.10.Jm, 65.40.-b, 72.15.-v 



I. INTRODUCTION 

Magnetic susceptibility x an d magnetization M mea- 
surements versus temperature T and applied magnetic 
field H have been widely used to determine the mag- 
netic properties of materials containing local magnetic 
moments. 1 Such measurements give valuable information 
about the magnetic interactions and magnetic structure 
of a material. In collinear antiferromagnetic (AF) ma- 
terials, the qualitative T dependence of the anisotropic 
X below the AF ordering (Neel) temperature Tn is well 
known. In 1941, Van Vleck calculated using molecu- 
lar field theory (MFT) the anisotropic x( T < t n) for 
a collinear AF containing identical crystallographically 
equivalent spins interacting by Heiscnbcrg exchange, but 
only for the special case of a two-sublattice bipartite 
spin lattice, i.e., in which the nearest-neighbor spins 
of a spin on one sublattice ("up spins") are members 
of the other sublattice ("down spins"). 2 He further as- 
sumed that spins on one sublattice only interact with 
the nearest-neighbor spins on the other sublattice and 
with the same strength. 2 Due to its limited applicability, 
few comparisons of experimental \{T < Tn) data with 
these theoretical predictions have been made. 

We recently formulated generic predictions using MFT 
of the anisotropic x(T < Tn) of both collinear and planar 
noncollincar AF structures for Heiscnbcrg spin systems 
containing identical crystallographically equivalent spins 
with arbitrary exchange interactions between arbitrary 
sets of spins. 3,4 Several comparisons of our theoretical 
predictions with experimental anisotropic \(T < Tn) lit- 
erature data for single crystals of known collinear and 
noncollincar AFs were made and reasonable agreement 



was found. 3 Such comparisons are expected to be most 
accurate for three-dimensional spin lattices with large 
spin S, which respectively minimize quantum fluctu- 
ations associated with low spin lattice dimensionality 
and/or low spin that are not taken into account by MFT. 
The MFT that we formulated is also expected to be 
most accurate for spin-only ions with angular momentum 
L = 0, which minimizes crystalline electric field effects 
arising from the spin-orbit interaction such as single-ion 
anisotropy effects. As discussed in Ref. 3, the deviation 
of the MFT prediction from the anisotropic %(T < Tn) 
data for Heisenberg spin systems can be used as a quanti- 
tative diagnostic for dynamical spin fluctuations and cor- 
relations beyond MFT. Comparisons of x(T) of an AF 
compound with MFT predictions have been used in the 
past to test for the occurrence of such dynamical short- 
range spin correlations, but usually only at temperatures 
above Tn- 

In Ref. 3 we pointed out that the same MFT pre- 
dictions used to extract information about the magnetic 
interactions and magnetic structures of AFs from anal- 
yses of anisotropic x(T < Tn) data for single crystal 
AFs should also be useful for analyzing the necessar- 
ily isotropic x(T < Tn) data for polycrystalline AFs. 
For example, such measurements can distinguish between 
collinear and planar noncollincar AF structures, even 
when multiple collinear AF domains occur. They can 
also be used to estimate the wave vector and turn an- 
gle between adjacent planes of spins along the helix or 
cycloid axis of planar helical or cycloidal AF structures. 5 

We report in this paper M(H, T) and x(T) mea- 
surements of polycrystalline EuNiGe3 and demonstrate 
that this compound exhibits long-range AF order below 
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T N = 13.6 K. We analyze the x(T) data at T < T N 
using our new MFT as well as at T > T N using the 
conventional Curie- Weiss law, which is also a MFT pre- 
diction. We also report x-ray diffraction measurements 
of the crystal structure of this material, electrical resis- 
tivity p{T) measurements fitted by the Bloch-Griineisen 
model and heat capacity C P (T) measurements analyzed 
using the Dcbyc model and correlate the results with the 
magnetic measurements. 

The compound EuNiGe3 crystallizes in the body- 
centered tetragonal BaNiSn3-type structure (space group 
74mm) with the Eu atoms at the corner and body-center 
positions forming a square lattice in the a6-plane stacked 
in a zigzag ABA fashion along the c-axis as shown be- 
low in Fig. 2(a). 6,7 No information is available about 
its physical properties. Recent measurements of the 
physical properties of the related compounds i?NiGe 3 
(R = Y, Ce-Nd, Sm, Gd-Lu) with the same stoichiom- 
ctry but with a different base-centered orthorhombic 
SmNiGe3-type crystal structure (space group Cmmm) 
have been reported, and most of these are found to or- 
der antiferromagnetically. 8 These compounds also con- 
tain a square R sublattice with a different ABBA zigzag 
stacking sequence, but which is qualitatively similar to 
the stacked Eu square lattice in EuNiGe3. Therefore 
EuNiGe3 also appeared to us to be a candidate for AF 
ordering as we subsequently confirmed. Also, spin-only 
Eu +2 ions with S = 7/2 have orbital angular momentum 
L = 0, an advantageous property resulting in negligible 
crystalline electric field effects and a spectroscopic split- 
ting factor g « 2. 

From analysis of our x{T) data for EuNiGc 3 at T > Tn 
using the Curie- Weiss law, the dominant interactions be- 
tween the Eu +2 spins S = 7/2 are ferromagnetic (FM), 
in spite of the collinear long-range AF order at T < Tn 
suggested by our x data at T < Tn. Taking into account 
the symmetry of the unit cell, we propose that the Eu 
spins within a tetragonal afc-plane interact ferromagncti- 
cally but spins in adjacent layers along the c-axis interact 
antiferromagnetically. We further propose that this set 
of exchange interactions gives rise to a collinear A-type 
AF structure in which the Eu ordered moments within a 
layer are aligned ferromagnetically with respect to each 
other, but are aligned antiferromagnetically with respect 
to the moments in the two adjacent planes along the c- 
axis. 

For the scheme of magnetic interactions in EuNiGe 3 
that we propose, the FM interactions between spins 
within an a6-planc act within the same sublattice, and 
hence the interactions are not consistent with Van Vleck's 
MFT for x(T < In) discussed above, even though both 
AF structures are collinear. This means that Van Vleck's 
prediction for x(T < Tn) is not appropriate for analyzing 
such data. Our generic MFT must be used instead. An 
analysis of x(T < ?n) for polycrystalline EuNiGe3 can- 
not determine the orientation of the easy axis of the A- 
type AF structure with respect to the crystal axes, e.g., 
along the c-axis or within the a&-plane, although such 



a determination is possible using x(T) data for single- 
crystal samples. Future magnetic neutron and/or mag- 
netic x-ray scattering measurements can test our model 
for the A-type AF structure and also determine the di- 
rection of the ordered moments. On the basis of the 
analysis of x(T > Tn) in terms of the Curie- Weiss law, 
we obtain estimates of the nearest-neighbor in-plane and 
out-of-plane Eu-Eu exchange interactions. 

The remainder of the paper is organized as follows. 
In Sec. II, we discuss the synthesis of the polycrystalline 
EuNiGe3 sample along with details of the measurements. 
The experimental results, analyses and discussion are 
presented in Sec. III. A summary and our conclusions 
are given in Sec. IV. 



II. EXPERIMENTAL DETAILS 

Polycrystalline EuNiGe3 was prepared from the high 
purity elements. Eu was obtained from Ames Labora- 
tory and Ni (99.996%) and Gc (99.9999+%) from Alfa 
Aesar. Eu pieces were surrounded by Ni and Ge pow- 
ders in a pressed pellet. The pellet was placed in a 2 mL 
alumina crucible and sealed in an evacuated fused silica 
tube. The sample was heated at 850 °C for 30 h followed 
by a thorough grinding to ensure homogeneity, and then 
after pelletizing heated in an evacuated fused silica tube 
at 900 °C for 7 d. All sample handling, except for brief 
exposures to air to press the pellets and load the silica 
tubes, was done in a glove box containing high-purity He 
gas. A single-phase sample (apart from a trace of elemen- 
tal Ge) was obtained as established from powder x-ray 
diffraction measurements described in Section III A be- 
low. Growths of single crystals using Sn and NiGe3 fluxes 
were attempted but were not successful. 

Powder x-ray diffraction (XRD) data were collected 
using a Rigaku Geigcrflcx diffractomctcr with Cu Ka ra- 
diation. Rietveld refinement of the XRD data was ac- 
complished using the FullProf package. 9 

M measurements versus H and T were carried out 
using a superconducting quantum interference device 
(SQUID) magnetometer (Quantum Design, Inc.). A gel 
cap was used as sample holder and its small diamagnctic 
magnetization was measured separately and corrected for 
in the magnetization data for EuNiGe 3 that are pre- 
sented. We use Gaussian cgs units for the magnetization, 
magnetic susceptibility and magnetic field throughout, 
where the magnetic field unit of Tesla, when it appears, 
is a unit of convenience (1 T = 10 000 Oe). 

Heat capacity C p and electrical resistivity p measure- 
ments were carried out using a Quantum Design Physical 
Property Measurement System (PPMS). The sample for 
C p measurements was attached and thermally coupled to 
the addenda with Apiezon N grease. The p measurements 
utilized a four-probe ac technique with the ac-transport 
option of the PPMS, where a rectangular parallelopipcd- 
shaped sample was cut from the sintered pellet for the 
measurements using a jeweler's saw. Platinum electri- 
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FIG. 1: (Color online) Room temperature powder XRD pat- 
tern (red circles) of EuNiGe3, Rietveld refinement fit (solid 
black line), difference profile (lower solid blue line), and posi- 
tions of Bragg peaks (vertical green bars). 




FIG. 2: (Color online) Comparison of the structure of (a) 
EuNiGe3 [BaNiSn3-type (space group 74mm)] with (b) the 
ThCr2Si2-type structure (space group 74/mmm). 



cal leads were attached to the sample using EPO-TEK 
P1011 silver epoxy and the sample was attached to the 
resistivity puck with GE 7031 varnish. The p(T) mea- 
surements were carried out on both cooling and heating 
to check for hysteresis. 



III. EXPERIMENTAL RESULTS, ANALYSES 
AND DISCUSSION 

A. Crystal Structure Determination 

The crystal structure of EuNiGc3 reported in Ref. 6 
was used as the starting point for the Rietveld refinement 



TABLE I: Crystal data for EuNiGe3 at room temperature 
(BaNiSn3-type structure: space group 74mm). Refined lat- 
tice parameters and unit cell volume are a = b = 4.3366(1) A, 
c = 9.8802(2) A and V cc ii = 185.81(1) A 3 . Listed are the 
Wyckoff atomic position and the atomic coordinates x, y, and 
z of each atom type. The goodness-of-fit parameters obtained 



are R p = 7.41%, 



10.4%, and X = 5.27. 
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Wyckoff position 
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Eu 


2a 








0.0028(3) 


Ni 


2 a 








0.6581(4) 


Gel 
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1/2 


0.2582(4) 


Ge2 


2 a 








0.4163(4) 



of our powder XRD data. We also observed very weak 
peaks from an impurity phase with the strongest peaks 
at diffraction angles 29 = 27.3° and 45.3°. These corre- 
spond to the strongest peaks of pure Ge, indicating the 
presence of a trace amount of elemental Ge in our sam- 
ple. During refinement of the XRD data for the EuNiGc3 
phase, the thermal parameters of the atoms were fixed at 
zero. In the final refinement the occupancies of the atoms 
were fixed at the stoichiometric values because no signif- 
icant difference in the goodness of fit was obtained when 
the occupancies were allowed to vary. 

A satisfactory Rietveld refinement of the powder XRD 
data for EuNiGe3 was obtained assuming the body- 
centered-tetragonal BaNiSn3-type structure (space group 
1 4mm) as shown in Fig. 1, with parameters listed in 
Table I. This structure and the common ThCi'2Si2- 
type structure are both derivatives of the BaAl 4 -type 
structure. 10 A comparison between the EuNiGc3 and 
ThCr2Si2 structures is shown in Fig. 2. The Eu and 
Th positions are identical in the two structures. The 
ThCr2Si2-type structure is centrosymmctric whereas the 
BaNiSn3-type structure is not. In the BaNiSn3-type 
structure, the transition metal square lattice in the ab- 
plane is rotated by 45° with respect to the a6-plane Cr 
square-lattice in ThCr 2 Si2, and the ordering of the Si or 
Ge layers and the transition metal layers along the c-axis 
is different. 

The refined lattice parameters for EuNiGe3 in the cap- 
tion of Table I can be compared with the reported values 
a = 4.737(2) A and c = 9.891(3) A. 6 Our a-axis parame- 
ter is much smaller by 0.400 A than the reported value, 
which we therefore assume is due to a typographical error 
in Ref. 6, but the c-axis parameters are nearly the same. 



B. Electrical Resistivity Measurements 

The p of EuNiGe3 was measured from T = 1.8 to 
350 K and the data are presented in Fig. 3. The sam- 
ple shows a very large residual resistivity ratio RRR = 
p(350 K)/p(1.8 K) = 54.8 indicating that our polycrys- 
tallinc sample is of good quality for electronic transport 
measurements. A strong decrease in p at T < Tn oc- 
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FIG. 3: (Color online) Electrical resistivity p versus temper- 
ature T for EuNiGe3 taken on heating (open red circles). For 
clarity, only every other data point is plotted. The black curve 
is a fit of the data at T > 14 K by Eq. (1). An extrapolation 
of the fit to T = is also shown. Inset: Expanded plot of the 
p(T) data at low T that were taken on cooling (black squares) 
and heating (open red circles). 



curs due to the loss of spin disorder scattering below 
Tn ~ 13.6 K (see also below), as shown on expanded 
scales in the inset of Fig. 3. The data in the inset also 
show no hysteresis between heating and cooling runs, in- 
dicating that the AF transition is thermodynamically of 
second order. Above Tn the resistivity due to spin dis- 
order scattering is expected to be constant. 11 Therefore, 
the T dependence above Tn is due to other electron scat- 
tering mechanisms. Typical mechanisms are electron- 
electron scattering which leads to a T 2 dependence and 
electron-phonon scattering with or without simultaneous 
Umklapp scattering. 

The Bloch-Griineisen model predicts the contribution 
to p(T) due to scattering of electrons by longitudinal lat- 
tice vibrations in the absence of Umklapp scattering. 12-14 
When additional constant terms are added to account for 
the residual resistivity (po) and the spin disorder resis- 
tivity at T > Tn (psd), the sum is 



p{T) = po + p s d 



T 
©R~ 



5 pBr/T 




(e x - 1)(1 - e- x ) 



^dx, (1) 



where 1Z is a material-dependent prefactor that is inde- 
pendent of T and Or is the Debye temperature deter- 
mined from resistivity measurements. 

To fit our p(T > T N ) data by Eq. (1), we uti- 
lized a high-accuracy analytic Pade approximant for the 
Bloch-Griineisen function in Eq. (1) that wc formulated 
recently. 14 As seen in Fig. 3, an excellent fit by Eq. (1) 
was obtained to the p(T) data for T > Tn- The parame- 



ters obtained from the fit are po + p s d = 10.21(3) pfl cm, 
6 R = 265(2) K and p(T = Or) = 44.3(2) /xficm. The 
quoted statistical errors on the resistivity contributions 
do not take into account an estimated systematic error 
of order 10% arising from uncertainty in the geometric 
factor and from the porosity and grain boundary scat- 
tering of the sintered sample. The excellent agreement 
of the temperature dependence of the data with the fit 
indicates that electron-phonon scattering is the primary 
scattering mechanism giving rise to the T dependence 
of p for T > Tn- This conclusion is supported by the 
agreement of Or = 265(2) K with the Debye tempera- 
ture 0d = 268(2) K obtained below by fitting the lattice 
heat capacity by the Debye model over approximately 
the same T range. 



C. Magnetization and Magnetic Susceptibility 
Measurements 

The magnetic susceptibility x = M/H versus T of 
EuNiGe3 was measured at H = 1 T and 3 T in the 
T range 1.8-300 K as shown in Fig. 4(a). The inverse 
susceptibility x^ 1 ^) for H = 3 T is plotted in Fig. 4(b), 
and x(T) data for H = 0.5 T from T = 1.8 to 25 K 
are shown in Fig. 4(c). The sharp peak at Tn ~ 14 K 
in the x( T ) data for if < 1 T in Figs. 4(a) and 4(c) is 
the signature of AF ordering occurring at this T. An- 
other transition at ps 5 K is observed as a sharp cusp 
for fields H < 500 Oe as shown in the inset of Fig. 4(a). 
The nature of this transition is unclear. There is no ev- 
idence of a transition at this T in either the p(T) data 
presented above or the C P (T) data presented below. A 
similar unexplained feature at about the same tempera- 
ture was observed in EuPdGe3. 15 The feature cannot be 
due to a ferromagnetic EuO impurity phase because the 
Curie temperature of EuO is 69 K. 16 

At T > Tn, the molar x follows the Curie- Weiss law 

x(t) = jSr p (2) 

as shown by the plot of x~ 1 versus T in Fig. 4(b) for 
H = 3 T, where C is the molar Curie constant and 9 P 
is the Weiss temperature. These x _1 (T) data decrease 
linearly with decreasing T down to w 50 K, below which 
short-range AF correlations begin to cause a deviation 
from the Curie- Weiss law. Thefore we fitted the data in 
Fig. 4(b) by Eq. (2) only from 50 to 300 K. The result- 
ing fit gave C = 7.67(2) cm 3 K/mol and p = 5.3(3) K. 
The C value is close to the value C = 7.88 cm 3 K/mol ex- 
pected for Eu +2 with S = 7/2 and g = 2. The agreement 
of the Curie constant with that expected for Eu +2 with 
S = 7/2 indicates that the Ni atoms are nonmagnetic. 
The conclusion that the Eu atoms have S = 7/2 agrees 
with the analysis of the magnetic entropy described in 
Section HID below. The ratio p /Tn is 
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FIG. 4: (Color online) Magnetic suceptibility \ ( a ) an d X _1 
(b) of polycrystalline EuNiGe3 versus temperature T. The 
inset in (a) shows expanded plots below 25 K of x(T) with 
H = 100 Oe to 3 T, and the straight black line in (b) is 
a fit of the data at T > 50 K by the Curie- Weiss law in 
Eq. (2). An extrapolation of the fit to the horizontal axis is 
also shown, (c) Expanded plot of x(T) at low T taken with 
H = 0.5 T. The curves are the MFT predictions for x±, Xll; 
and Xpowder versus T from Eqs. (8) and (9) compared with 
the experimental Xpowder data (open red circles). 



where we used the precise value Tn = 13.6 K determined 
from our heat capacity measurements below. 

Using the Heisenberg Hamiltonian H = JijSi-Sj 
for a system comprised of identical crystallographically 
equivalent spins, where the sum is over distinct pairs of 
spins, from MFT one can write 9 P and Tn in terms of the 
exchange interactions Jij between spin i and its neigh- 
bors j as 3 



S(S + l) 
3fc B 

S(S + 1) 
3fc B 



3 



cos < 



(4a) 
(4b) 



where /cb is Boltzmann's constant and <pji is the angle 
between ordered moments j and i in the magnetically- 
ordered state. There is no restriction on the range of the 
exchange interactions Jy in Eqs. (4), and these can there- 
fore be nearest-neighbor, next-ncarcst-ncighbor, etc., in- 
teractions. 

From Eq. (4a), the positive value of 9 P observed for 
EuNiGc3 indicates that the dominant interactions be- 
tween the Eu spins are FM (negative), in spite of the long- 
range AF ordering. In order to simultaneously satisfy 
these two conditions within the symmetry constraints of 
the crystal structure, we propose that the dominant FM 
interactions J\ are between nearest-neighbor Eu spins 
in the a6-plane, with subdominant AF interactions J c 
between nearest-neighbor Eu spins in adjacent layers. 
These interactions give rise to an often-observed A-type 
AF structure in which FM aligned layers of ordered Eu 
moments in the a6-plane are AF aligned with the Eu 
moments in adjacent Eu planes along the c-axis. Our 
magnetization data for polycrystalline EuNiGe3 cannot 
determine the axis along which the ordered moments arc 
aligned, but single-crystal measurements could determine 
that. A-type AF structures were reported for the mag- 
netic Co atoms in single-crystal CaCo2As2 with the Co 
ordered moments aligned along the c-axis, 17-19 and for 
the magnetic Eu atoms in single-crystal EUFC2AS2 with 
the Eu ordered moments aligned in the et6-plane. 20 Both 
compounds have the tetragonal ThCr2Si2-type crystal 
structure shown in Fig. 2(b). 

The Weiss and Ncel temperatures for an A-type AF of 
stacked square lattices as in Fig. 2(a) with only nearest- 
neighbor interactions are given by Eqs. (4) as 

9 p = - ^fr-^ ( 4J i+ 8J e), (A type AF) (5a) 



3/cb 
S(S+1) 
3 kvt 



(4Ji - 8J C ), 



(5b) 



where for the A-type stacked square lattice AF, by def- 
inition one has <pji = for the four nearest-neighbor 
spin pairs within an ab plane and <pji = 180° for the 
eight nearest-neighbor spin pairs between adjacent lay- 
ers along the c-axis. In the latter case, due to the lack 
of a horizontal mirror plane through the body-centered 
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Eu site in the noncentrosymmctric crystal structure in 
Fig. 2(a), the two J c values from an Eu spin to the four 
nearest-neighbor Eu spins in each of the two adjacent ab- 
plane layers, respectively, are different. Therefore the de- 
rived J c is an average of the two interplanar interactions. 
For a spin lattice consisting of square lattices stacked di- 
rectly above and below each other, the coefficient of J c in 
Eqs. (5) would have been 2 instead of 8. From Eqs. (5) 
one can solve for the two exchange interactions J\ and J c 
in terms of the measured values of Tn and 9„, yielding 



3(T N + p ) 



85(5 + 1) 
J c 3(Tn — P ) 
Afe ~ 165(5 + 1)' 



, (A type AF) (6a) 
(6b) 



Using T N = 13.6 K, 6 P = 5.3(3) K and 5 = 7/2, Eqs. (6) 
yield 

= -0.45(1) K, = 0.099(4) K. (7) 

Kb kb 

These results quantitatively confirm the above qualita- 
tive deduction based on the positive Weiss temperature 
that the dominant Eu-Eu exchange interactions in the 
system are ferromagnetic (negative). These we deduce 
to be the J\ interactions between Eu spins within an 
a6-planc layer, whereas the inter layer interactions J c are 
antifcrromagnetic (positive). In particular, the dominant 
Ji (i = 1 or c) value as well as the dominant Zj Ji value are 
both negative, where Zi = 4 and 8 are the coordination 
numbers of Eu by Eu for in-plane J\ and out-of-plane J c , 
respectively. 

For a Heisenberg system of identical crystallograph- 
ically equivalent spins in the absence of magnctocrys- 
tallinc anisotropy, our MFT predicts the anisotropic tem- 
perature dependence of the susceptibility at T < Tn for 
a collinear AF such as the A-type AF to be 3 



X±(T < T N ) 

X|| (*) 
(5 + 1)* 



1-/ 
*(*)-/ 



yo 



3/i 



(S + l)t' 



(8b) 
(8c) 



with the Brillouin function Bs(y) and its derivative 
B's(y) respectively given by 



Bs(y) 
B' s (y) 



- 

= dB s (y) 
dy 



l)coth (25 



l) V - 

'2 



-coth(f)}, 
(8d) 



csch 2 (;y/2) - (25 + l) 2 csch 2 [(25 + l)y/2] 
45 : 



where || and _L refer to the magnetic field applied paral- 
lel and perpendicular to the easy axis, respectively, the 
reduced temperature is t = T/T$, and we use the un- 
conventional definition of B${y) in Rcfs. 4 and 21. The 



reduced T-depcndcnt ordered moment is pQ = /io/^sat 
where [1$ is the magnitude of the ordered moment at 
H = and ^ sa t = gS/iB is the saturation moment. 
The fto(t) is determined by numerically solving the ex- 
pression /2o = Bs(yo)- 3,4 Equations (8) predict that 
X\\(T = 0) = and that the susceptibility is isotropic 
at T = T N , i.e., X||(T N ) = X± (T N ) = x(T N ). 3 The X fol- 
lows the Curie- Weiss law in Eq. (2) and is isotropic for 
T > T N . 3 < 4 

In a polycrystalline sample such as ours, it is assumed 
that the many small crystallites are randomly oriented. 
Therefore, the x(T) can be obtained as the spherical 
("powder") average of the x± an d X|| components to be 



X P owdci-{T) = -[2x_l(T) +%||(T)] 



(9) 



This powder average is the same if multiple A-type AF 
domains occur. Using the values 5 = 7/2, / = 0.39 from 
Eq. (3) and the observed x(^n) = 0.66(1) cm 3 /mol from 
Fig. 4(c), the T dependence of the powder average sus- 
ceptibility for T < Tn obtained from Eqs. (8) and (9) is 
the green line in Fig. 4(c) with no adjustable parameters. 
The predicted T dependence for T < Tn in Fig. 4(c) is 
in reasonable agreement with the measured x(T) data 
(open red circles) plotted in the same figure. The exper- 
imental x(T — > 0) value in in Fig. 4(c) is slightly lower 
than the MFT prediction, which may arise from a slight 
deviation from a random distribution of the orientation 
of the grains in the polycrystalline sample. 

M versus H isotherms were measured for EuNiGe3 in 
the H range 0-5.5 T at various temperatures as shown in 
Fig. 5. The MFT prediction per spin in the paramagnetic 

r N 



Hz 



.(5 + 1)* 



(10) 



where ji z is the magnetization induced in the direction 
of the applied field by the applied field, fl z = fi z / ^ sat 
and the reduced applied magnetic field h is defined as 
h = g^H I '(/cbTn). AH of the parameters g, 5, / and Tn 
in Eq. (10) were already determined above. The M(H) 
isotherms for T > Tn calculated from numerical solution 
of Eq. (10), where M = N^, z and N is the number of 
spins, are compared with no adjustable parameters with 
the corresponding experimental M{H) data in Fig. 5. A 
proportional M(H) behavior is predicted and observed 
for T > 50 K, whereas negative curvature in M{H) is 
predicted and observed at 20 K. However, the calculated 
curve for T = 20 K is slightly above the observed data 
due to dynamical short-range AF ordering in the sample 
on approaching Tn = 13.6 K from above, which sup- 
presses the magnetization. Such dynamical short-range 
ordering effects above Tn are not taken into account in 
MFT (see also the next section) . At temperatures below 
Tn, the maximum observed magnetization of 6.70 /ie/Eu 
at T = 1.8 K and H = 5.5 T is approaching the satu- 
ration moment /i sat = 7 ^e/f.u. expected for Eu +2 with 
5 = 7/2 and g = 2. 
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FIG. 5: (Color online) Magnetization M versus applied field 
H isotherms measured at the indicated temperatures for 
EuNiGe3. Solid curves of corresponding color are the the- 
oretical predictions by MFT in Eq. (10) for the data in the 
paramagnetic state (T > Tn — 13.6 K). In the ordinate label, 
f.u. means formula unit. The deviation of the fit from the data 
for T = 20 K is due to the presence of short-range AF order 
at that T, and the S-shaped behavior for T < Tn is due to 
occurrence of spin-flop transitions with a distribution of spin 
flop fields due to the poly crystalline nature of the sample. 



The M(H) isotherms below T N at T = 1.8, 5 and 10 K 
show an S-shaped dependence on H. Qualitatively, this 
can be explained by a scries of field-induced first-order 
spin-flop transitions where the ordered moments flop to 
a perpendicular orientation with respect to the applied 
field. In order for a spin flop transition to occur in a 
collinear AF, some type of magnetocrystalline anisotropy 
must be present that aligns the moments along the easy 
axis at zero field. The first-order spin flop transition docs 
not occur at a single field as observed in a single crystal 
with the field along the easy axis because of the ran- 
dom orientations of the crystallites in the polycrystallinc 
sample. In a polycrystallinc sample, one expects the spin 
flop field in a grain with its easy axis at an angle 9 to 
the field to obey H Rop (6) = H &op {6 = 0)/cos<9. Thus 
the spin flop field increases with increasing 9. When the 
easy axis is perpendicular to the applied field (9 = 90°), 
a spin flop transition is not possible because the orien- 
tation of the ordered moments is already perpendicular 
(in H = 0) to the field direction. From Fig. 5, we infer 
that Fflop(6> = 0, T 0) ~ 1.5 T. A calculation within 
MFT of the powder- averaged M(H) for a polycrystallinc 
sample, incorporating both the anisotropy field and the 
^-dependent distribution of spin flop transition fields, is 
beyond the scope of the present work. 
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FIG. 6: (Color online) Heat capacity C p versus temperature 
T for EuNiGe3 . The red curve is a fit by a Pade approximant 
of the Debye function in Eq. (11) to the data at T > 25 K. 
An extrapolation of the fit to T = is also shown. Inset: 
Expanded plot of C P (T) for EuNiGe 3 at low T (black cir- 
cles), together with background C P (T) data for LaNiGe3 (red 
squares) 22 with a renormalized T scale (see text). 



D. Heat capacity measurements 

The C p of EuNiGe3 was measured at H = in the 
T range 1.8-300 K and the data are plotted in Fig. 6. 
A sharp A-shaped peak is observed at Tn = 13.6 K as 
shown in more detail in the inset of Fig. 6, confirming 
that the AF transition observed in the magnetization 
measurements is a bulk magnetic phase transition. The 
C p (300 K) = 122 J/(molK) is approaching the classical 
Dulong-Petit high-T limit Cy = 3nR = 124.7 J/(mol K) 
for the heat capacity of acoustic lattice vibrations at con- 
stant volume, where n = 5 is the number of atoms per 
formula unit and R is the molar gas constant. 

The Debye model describes the heat capacity versus T 
due to such lattice vibrations by 1 



I 7T- 

fcln 



3 ,e D /T 



x 4 e x 



■ dx, 



(11) 



where 0d is the Debye temperature determined from 
heat capacity measurements. In addition, for a metal 
one can add a linear jT term to Eq. (11) to account for 
the electronic specific heat contribution, where 7 is the 
Sommerfeld electronic specific heat coefficient, and for a 
magnetic material one can add the magnetic contribution 

C m ag(T). 

An accurate analytic Pade approximant of the Debye 
function that we recently formulated to simplify fitting 
of experimental C p (T) data by the Debye theory 14 was 
used in place of Eq. (11) to fit our data. We fitted our 
C P (T) data from T = 25 to 300 K because below sa 25 K 
the magnetic heat capacity contribution C mag becomes 
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significant (see below). Also, because of the presence of 
C mag , 7 could not be accurately determined from the 
C P (T — ► 0) data. When allowed to vary, it refined to 
the value 2(2) mJ/molK 2 . Therefore, 7 was fixed at zero 
for the final fit. The only adjustable parameter in the 
final fit was 9d, which was found to be Od = 268(2) K. 
As seen in Fig. 6, a reasonably good fit of the data by 
the Dcbye model is obtained over the entire temperature 
range above 25 K. 

In order to isolate C mag (T), C p (T) data for the 
isostructural nonmagnetic reference compound LaNiGe 3 
from Ref. 22 was used as a heat capacity background, af- 
ter correcting for the formula weight (FW) difference by 
multiplying the temperature scale for C p (T) of LaNiGc3 

by y/FW LaNiGc3 /FW E uNiGo 3 ■ Fr om the inset of Fig. 6, 
the renormalized C P (T) for LaNiGe 3 is seen to be a 
reasonable estimate of the background heat capacity of 
EuNiGe3. The C mag (T) obtained by subtracting the 
renormalized C P (T) of LaNiGe3 from the C P (T) data for 
EuNiGc 3 is plotted in Fig. 7(a) and C mag (T)/T is plotted 
in Fig. 7(b). 

MFT predicts that C mag of a Heisenberg spin system 
containing identical crystallographically equivalent spins, 
per mole of spins, is 4 

ir = ~sTT M<) ~tf~- (12) 

There are no adjustable parameters in this prediction if S 
and Ik are known, as they are here. Comparisons of the 
predicted C mag (T) and C mag (T)/T from Eq. (12) with 
our experimental data are shown in Figs. 7(a) and 7(b), 
respectively. The hump in the experimental C mag (T) at 
T ~ In/3, which is much more pronounced in the plot of 
C ma g(T)/T 1 , is reproduced by the MFT calculation. This 
hump in C mag (T) increases in magnitude as S increases, 
is particularly noticablc for S = 7/2, and is not visible 
for, e.g., S = 1/2. 4 The hump arises in MFT from the 
combined effects of the T dependence of the energy split- 
ting of the Zeeman levels of the spin due to the T depen- 
dence of the ordered moment and resulting T dependence 
of the exchange field seen by each moment, together with 
the T dependence of the Boltzmann populations of those 
levels. The hump in C mag (T) observed at T ~ Ik/3 
for magnetically ordered compounds containing Eu +2 or 
Gd +3 with S = 7/2 is sometimes misinterpreted as aris- 
ing from either an electronic magnetic Schottky anomaly 
that is combined with a T 3 magnon contribution to re- 
produce the observed C mag (T), as evidence for some type 
of magnetic phase transition, or as a giant nuclear Schot- 
tky anomaly induced by the ordered moments. 

The magnetic contribution S mag (T) to the entropy was 
calculated from the C mag (T) derived from our experi- 
ments according to 

S mas (T) = £ dT. (13) 

Because C mag (T) data were not obtained below 1.8 K, 
the C mag /T data were extrapolated from T = 1.8 K to 
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FIG. 7: (Color online) (a) Magnetic contribution C mag to the 
heat capacity of EuNiGe3 versus temperature T. (b) C mag /T 
versus T. The black curve at low T in (b) is an extrapolation 
of the data from T = 1.8 K to T = 0. The red curves in 
(a) and (b) are the predictions of MFT in Eq. (12) for spin 
S = 7/2. (c) Magnetic contribution S mag to the entropy 
versus T. The horizontal red line is the expected high-T limit 
SWag = i?ln(8) = 17.29 J/molK for S = 7/2. 



T = using, for simplicity, the T 2 dependence predicted 
by spin wave theory for a three-dimensional AF in the 
absence of an anisotropy gap. The calculated entropy 
between T = and 1.8 K on the basis of this extrap- 
olation is 0.30 J/mol K. This is an upper limit since 
the presence of an anisotropy gap would instead give 
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an exponential decrease in C mag (T) below 1.8 K. As 
seen in Fig. 7(c), the molar S mag saturates to the value 
Smag(T — > oo) ps i?ln(8) expected from quantum statis- 
tics according to S l mag (T — > oo) = R\n(2S + 1), where 
S = 7/2 for Eu+ 2 . Even though the MFT prediction has 
significant deviations from the measured C mag /T data at 
T < Tn in Fig. 7(b), the area (magnetic entropy) between 
the calculated and observed data that is missing below 
T N is recovered above T N . The finite C mag (T) for T > T N 
arises from dynamical spin correlations at T > Tn , as of- 
ten observed, that are not taken into account by MFT. 
As a result, S niag (T N = 13.6 K) = 14.3 J/mol K is 83% 
of the asymptotic high-T limit, with the remainder re- 
covered by « 30 K = 2.2 T N . 



IV. SUMMARY AND CONCLUSIONS 

A nearly single-phase polycrystalline sample of 
EuNiGe3 was synthesized and its physical properties were 
investigated. Rietveld refinements of the powder XRD 
data confirmed that this compound crystallizes in the 
body-centered-tetragonal BaNiSn 3 -type structure with 
space group IAmm as previously reported. The p, x 
and C p measurements consistently reveal an AF ordering 
transition at Tn = 13.6 K. A summary of some of the 
results from these measurements is given in Table II. 

The p(T) measurements of EuNiGe3 reveal metallic be- 
havior. The large RRR = 54.8 indicates the high quality 
of the sample. The data for T > Tn are well-described 
by the Bloch-Grunciscn theory for the T-dependent re- 
sistivity arising from electron-phonon scattering. A fit 
to the data for T N < T < 350 K by the theory yielded 
a Debyc temperature Or = 265(2) K. The p decreases 
rapidly on cooling below Tn due to loss of spin-disorder 
scattering. Since the compound is metallic, the magnetic 
coupling between the Eu spins likely arises mainly from 
the indirect RKKY interaction mediated by the conduc- 
tion electrons. 

The M(H,T) and X {T) measurements of EuNiGe 3 
showed the presence of long-range AF order in this sys- 
tem at T N wl4K. A fit of x -1 ^) by the Curie-Weiss 
law at T > 50 K revealed a Curie constant consistent 
with the presence of Eu +2 ions with S = 7/2 and .9 = 2, 
and a positive Weiss temperature 9 P = 5.3(3) K, indi- 
cating that ferromagnetic interactions are dominant de- 
spite the occurrence of long-range AF ordering. There 
was no evidence from our measurements that the Ni 
atoms are magnetic. A low-field xC^O measurement at 
T < Tn was compared with our prediction 3 from MFT 
for the polycrystalline average of the anisotropic x(T) of 
a collinear antiferromagnet below its Neel temperature 
and good agreement was found. Carrying out such a fit 
for a polycrystalline AF was one of the goals of this work 
as discussed in the introduction. A field-induced spin flop 
transition was inferred from the S-shaped M(H) curves 
at T < T N , with a low-T onset field of H ~ 1.5 T that 
was spread out to higher fields due to the polycrystalline 



nature of the sample. 

The presence of dominant FM interactions in EuNiGe3 
that orders antifcrromagnctically led us to propose that 
these interactions are between Eu spins within the ab- 
plane, with subdominant AF interactions between spins 
in adjacent planes along the c-axis. From these inter- 
actions, we propose that the collinear AF structure is 
A-type, in which the Eu spins within an a&-plane in 
Fig. 2(a) are aligned ferromagnetically with respect to 
each other and the spins in adjacent layers along the c- 
axis are aligned antifcrromagnctically to each other. The 
most likely ordered moment axis is cither the c-axis or 
an axis in the a&-plane. If it is the c-axis, magnetocrys- 
talline anisotropy effects would presumably not cause a 
distortion of the crystal structure on cooling below Tn- 
However, if the ordered moments are in the a6-plane, 
the A-type collinear ordering breaks the fourfold rota- 
tional symmetry about the c axis of the tetragonal room- 
temperature crystal structure, and an orthorhombic crys- 
tal distortion may be expected to occur on cooling below 

Our C P (T) data were fitted by the Debye model from 
T = 25 K to 300 K, yielding a Debyc temperature Od = 
268(2) K. This value is the same within the error bars as 
the Debye temperature Or = 265(2) K determined from 
the p(T) measurements, a rare occurrence. 14 The data 
exhibited a sharp A-shaped peak at Tn, which allowed 
the precise value of the Neel temperature to be deter- 
mined to be Tn = 13.6 K. The magnetic heat capacity 
contribution C mag (T) and the magnetic entropy S ma ^(T) 
were extracted and analyzed by MFT. The high-T lim- 
iting value of the entropy S mag ps i?ln(8) = Rhi(2S + 1) 
is consistent with our \{T > ^n) data that indicated 
S = 7/2. Significant short-range AF correlations occur 
above Tn, with about 83% of the maximum magnetic en- 
tropy present at Tn and the remaining 17% recovered by 
w 30 K = 2.2 T N . 

As shown in Fig. 2, EuNiGe3 crystallizes in the 
BaNiSn3-type structure which is similar to the ThCi - 2Si2- 
type structure. In fact, the compound EUM2GC2 crys- 
tallizes in the latter structure and its properties have 
been measured. 23,24 The Eu sublattices in the two com- 
pounds arc identical and the Ni atoms in both com- 
pounds are believed to be nonmagnetic. Therefore, we 
compare some properties of these two compounds in Ta- 
ble II. Both compounds contain Eu +2 ions with S = 7/2 
and g = 2 and order antiferromagnetically with EuNiGe3 
having the lower Tn- A significant difference between 
these compounds is the positive 6 p in EuNiGc3, indicat- 
ing dominant FM interactions as discussed above, and 
a negative one in EuNi2Ge2 indicating dominant AF in- 
teractions. This difference indicates that the magnetic 
interactions between the Eu spins are quite different in 
the two compounds and therefore suggests that the resul- 
tant AF structures may also be different. As discussed 
above, we propose that EuNiGe3 has an A-type AF struc- 
ture. The anisotropic xC^O measurements on EUM2GC2 
single crystals 23 suggest that the ordered moments lie 
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TABLE II: Summary of the physical properties of EuNiGe3. The properties of EuNi2Ge2 are also presented for comparison. 
Listed are: the tetragonal lattice parameters a and c at room temperature, the Neel temperature Tn, Weiss temperature 8 P , 
Curie constant C, effective magnetic moment fx e g = v8C of the Eu, and Debye temperature determined from heat capacity 
(0d) and resistivity (0r) measurements. 



Compound 


a 


c 


T N 


9 P 


C 


HeS 


e D 


6r 


Ref. 




(A) 


(A) 


(K) 


(K) 


(cm 3 K/mol) 


(Mb) 


(K) 


(K) 




EuNiGe 3 


4.3366(1) 


9.8802(2) 


13.6 


5.3(3) 


7.66(2) 


7.83(2) 


268(2) 


265(2) 


This work 


EuNi 2 Ge 2 






30.8 


-9.7 




7.69 






23 


EuNi 2 Ge 2 


4.144(3) 


10.15(1) 


30 


-8 


7.7 








24 



in the afr-plane with a collinear AF structure and multi- 
ple AF domains. Alternatively, a comparison of our re- 
cent predictions 3 of the anisotropic x(T) of planar non- 
collinear AFs with the x{T) data 23 for EuNi 2 Gc2 sug- 
gests that this compound may have a planar noncollinear 
AF structure with the ordered moments aligned within 
the a&-plane. It would be useful and interesting to deter- 
mine the AF structures of both EuNiGe3 and EuNi 2 Ge 2 
by magnetic neutron or x-ray scattering measurements 
and to correlate the results with the respective x(T) data 
for these compounds. 
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